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ABSTRACT 



A new model of QCD parton shower is proposed which is dedicated to two-photon 
process in e + e~ scattering. When hadron jets are produced, the photon may resolve 
into quark- ant iquark pairs so that the structure functions of the photon should be 
introduced. Based on the Altarelli-Parisi equation for these functions, an algorithm 
is formulated that allows us to construct a model for parton showers for the photon. 
Our model consists of two parts, one of which describes the deep inelastic scattering 
of the photon and the other one the scattering of two quasi-real photons. Using the 
model some results are presented on parton distributions and jet production. 
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Section 1. Introduction 

Recent extensive studies in TRISTAN experiments on the two-photon process in 
e + e~ scattering, particularly on jet production |T], have stimulated more detailed 
theoretical investigations in QCDf||, f|, [5|. However these calculations have been 
limited to the inclusive case so far. In order to make a precise comparison between 
experiments and theory, one needs an event generator where one starts with the 
generation of partons and then convert them into hadrons and other stable particles. 
Distributions of generated partons are completely determined by perturbative QCD. 
To generate partons one has two possibilities. The first one is the matrix element 
method which uses the fixed order perturbative calculation and the other one is the 
parton shower method|| based on the Altarelli-Parisi(AP) equation. In the latter 
one can take into account all the leading contributions from the collinear singularities 
without limitation in the order of the coupling. 

The parton shower model has been demonstrated to be powerful and indispens- 
able for quantitative studies in e + e~ annihilation and other processes. It is also true 
for the two-photon experiments and it is highly desirable to develop an event gener- 
ator that includes parton showers, since recent experiments provide data with rather 
good accuracy. In this paper we construct such a parton shower model dedicated to 
the two-photon process. 

One of the main features of this parton shower is the resolution of the photon, 
i.e. the case where the photon splits into a quark- ant iquark pair at the beginning 
of the reaction. Partons inside the photon have two sources. One is taken into 
account by the so called vector meson dominance model(VDM) and the other one 
is the perturbative photon. They correspond to the solutions of the homogeneous 
and inhomogeneous AP equations, respectively. In terms of parton language, the 
VDM contribution arises when the photon converts into the vector mesons like p° 
and u and then the partons inside these mesons undergo QCD evolution. On the 
other hand the perturbative photons are those that decay into quark-antiquark pairs 
and then these quarks evolve according to perturbative QCD. We will call them the 
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VDM and the perturbative photon parts in the following. 

In the case of two quasi-real photon scattering we have another contribution 
where the original photon undergoes a hard scattering. This shall be called the 
direct part. 

The present model is limited to the leading-logarithmic(LL) approximation. An 
extension to the next-to-leading logarithmic case is, however, not difficult 0. One 
reason for this limitation is that experimental data are not yet precise enough to 
determine the QCD scale parameter A and another reason is that the most important 
contributions of the higher order are already included into the model as discussed 
in0. 

To apply the model to actual reactions one has to distinguish two situations(see 
Fig.l). The first one is the deep inelastic scattering of a photon, i.e. one photon is 
quasi- real while the other has a large virtual mass squared — Q 2 . One requirement 
from perturbative QCD is that this virtual mass squared is much larger than the 
QCD scale parameter and hadron mass scale. Hence one has to assume Q 2 > 
5-10 GeV 2 . The second case is for two quasi-real photon scattering. In this 
process one can get reliable predictions from QCD calculations only when jets with 
large p\ are produced. Notice that our event generator produces only partons, not 
hadrons. Hence the VDM contributions express only the part related to the hard 
scattering. We do not consider the low pt scattering in the hadron interaction. 

In the next section we show how to formulate the parton shower model, starting 
from the AP equations. In section 3 a detailed description of the model will be 
presented. There we discuss the differences between the models of deep inelastic 
scattering and two quasi-real photon scattering. These are related not to the algo- 
rithm of the model but the way to apply it. Section 4 shows some results obtained 
by Monte Carlo simulations of the parton shower. One will see that it reproduces the 
same distributions as the conventional QCD calculations for inclusive observables. 
The last section is devoted to summary and discussion. 
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Section 2. Basic formulation 

The algorithm of the parton shower is based on the AP equation, which governs 
the evolution of the photon structure function. First we solve this equation in the 
way most convenient to our purpose, which gives us a clear basis for the construction 
of the shower model. For simplicity we limit the discussion here to the non-singlet 
case within the LL approximation. 

The AP equation for qNs(%, Q 2 ), the non-singlet quark distribution in the photon|| 
U, is given by 

Y2\ 



dq NS (x,Q 2 ) _ a s ri dy (0) 2 a 

d\nQ* -*kL 7 qq {x/y)qNs{y ' Q > + 2^ 



(0) (x) 



(1) 



Here P$ (x) is the quark splitting function. This equation differs from the usual 
one by the existence of an inhomogeneous term, k^ s (x). It is known that this term 
has a x-dependence of the form x 2 + (1 — x) 2 . The QED coupling is denoted as a 
and that of QCD is given by 

47r _ a 2 



a s (Q 2 



A) = n 



-N, 



(2) 



A)hiQ 2 /A 2 lnQ 2 /A 2 ' 3 
where Np is the number of flavors. Eq.(l) can be solved easily if one takes its 
moments: 

qNs(n,Q 2 )=[ dxx^qNsix.Q 2 ), (3) 
Jo 
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where d(n) is the moment of P^(x)/2n and a NS (n) is that of (a/2n)k 



(o) 

NS 



X) 



In order to show the dependence on n in a more transparent form, we introduce 
an integral to rewrite the second term, 

q N s(n, Q 2 ) = q N s(n, Q 2 )e a ° d ^ + -^—a NS {n) f dr,e^e a ^\ (5) 

®s{Q ) Jo 
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where 

3 = ln ««M) =ln M0W 

1,1 <*,(<« ln ln(Q§/A2)' W 
This is the final form of the solution in the moment space. 

To get the corresponding solution in the x-space, we make the inverse Mellin 
transformation defined by the formula 

/(*) = / ^~ n f(n), (7) 

where tq is a real number which defines the location of the integration path. The 
result can be cast into the following form 

Qns(x,Q 2 ) = f —K§l(x/y,s)q NS {y,Ql) 
Jx y 

+ \n(Q 2 /A 2 ) ['dT,e-* t ^Kf^x/y^ kf s {y). (8) 
Jo Jx y 2tc 

Here we have introduced the QCD kernel function K^ s (x, s) defined in ||. 

fn\ rro+ioo A n 

K$l(x/y,8)= ™-n e *o<i(n)s (g) 

Jr^—ioo Z7TI 

This kernel function is parametrized in a simple form so that it is easy to carry out 
convolution integrals in Eq.(|]) numerically. The physics meaning of each term in 
this solution is clear. The first term expresses the contribution from VDM, while 
the second contains the effect of the perturbative photon. 

In order to translate the Eq.(H) into the language of parton showers, it is conve- 
nient to change the integral variable r\ as follows, 

, lnQ 2 /A 2 , s 

where — K 2 can be regarded as the squared virtual mass of the quark that is emitted 
from the photon. By this change of variable, one can rewrite the solution in the 
following way: 

- 1 dy (0) 



qNs(x,Q 2 ) = I —KP s (x/y,s)q NS (y,Ql) 
Jx y 



Q'dK 2 [idy rX o), , 



The second term of this equation says that at the beginning the photon resolves into 
a quark-antiquark pair and then the spacelike quark(antiquark) evolves according to 
the usual QCD branching, because K^ s (x/y,n(K 2 )) expresses the evolution from 
K 2 to Q 2 . Based on this observation we can formulate an algorithm for the QCD 
evolution of the photon. So far the discussion is limited to the non-singlet case, but 
the extension to the singlet case can be done in a parallel way and our model is, 
needless to say, constructed to contain both contributions. 

Section 3. The model 

In this section we describe the model in some details. We assume the forward 
evolution method for the spacelike QCD evolution, as there is no essential difference 
between the forward evolution and the backward one p. One reason is, however, 
that the former makes it easier to write programs according to Eq.(|T|). 

The photon branching differs slightly from the other processes such as e + e~ an- 
nihilation or deep inelastic ep scattering. In the latter the total energy is conserved 
and the anomalous dimension corresponding to it vanishes identically. In the shower 
model which has been developed by one of the authors, the algorithm crucially re- 
lies on this conservation |n|. In other words the evolution is considered not for the 
particle distributions q(x,Q 2 ), but for the energy distributions xq(x,Q 2 ). In the 
present case, however, Eq.(|l|) is inhomogeneous and the second moment of k^ s (y) 
does not vanish so that the energy sum rule is broken. Hence in the photon shower 
model, one has to abandon parton energy conservation. Of course the total energy 
including that of the photon should be conserved. 



Therefore the quantity appearing on the right-hand side of this equation is taken 
into account as the weight of an event. 



jf dxx[^(q f (x, Q 2 ) + q f (x, Q 2 )) + G(x, Q 2 )} 




(12) 
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Let us state the algorithm step by step. It can be divided into the spacelike 
evolution and the hard scatterings that take place afterwards. Procedures of the 
spacelike evolution can be summarized as follows; 

(A) Choose Q 2 . 

(B) Calculate the energy in the VDM, which is independent of Q 2 , and the energy 
in the perturbative photon by using Eq. (|T2D . The sum of these energies is kept as 
the weight of an event. 

(C) Determine the process, either the VDM or the perturbative photon, according 
to the ratio of energies. 

(D) If the VDM is chosen, make the usual QCD evolution up to Q 2 from Qq, the 
cut-off momentum. 

(E) In the case of the perturbative photon, determine the virtual mass squared 
K 2 according to the probability dK 2 /K 2 . Then choose a flavor of quark according 
to the ratio of charges squared. 

(F) For the selected flavor of quark or antiquark, make the usual QCD evolution 
up to Q 2 from K 2 . 

The same algorithm is used for the deep inelastic photon scattering and two 
quasi- real photon scattering with large p 2 -,. The hard scattering parts, however, 
completely differ from each other. 

First we discuss the deep inelastic case. Here the QCD evolution gives the struc- 
ture function F 2 (x, Q 2 ) of the photon, which is connected to the cross section in e-7 
scattering through the following relation H, 

da Ana 2 ,„ . _, . 

y = Q 2 /(xs 7 ), 

where s 7 is the energy squared of the e-7 system. 

In the hard scattering of this process the gluon does not participate but the quarks 
do with the weight of the squared charge. Determination of the four momenta is 
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completed after the timelike evolution is developed and the energy of the initial 
photon is fixed by the spectrum function of the photon in the electron, which is 
given by the Weizsacker- Williams approximation. One of its simplest form is given 
by 



Fy/e(z,E) = -(l + (l- ^) 2 )(ln— -i). 

txz m e 2 

In the analysis of actual data more sophisticated versions are used. 

Next we discuss the scattering of two quasi-real photons. In this case contri- 
butions from the direct photon should be taken into account. The magnitude of 
the cross sections for various hard processes is necessary to determine the weight of 
events. They are different depending on the type of partons involved in the hard 
scattering. For example, if the initial state of the hard scattering is a gluon and a 
direct photon, the final state is a quark and an antiquark. Then we have to know 
the cross section. 

^(G + 7-, + S ). d3) 

Here s is the energy squared of this hard scattering and the hard scale Q 2 can be 
identified as the transverse momentum squared of the scattered partons, Q 2 = r p\. 

As the possible initial states of the hard scatterings we have to consider q-q, q-q, 
q-G, G-G, q-j, G-j and 7-7. 

In order to generate events in the most effective way with 0(1) weights, we adopt 
the following method. The cross sections for hard scattering with photons of direct, 
direct-resolved, doubly- resolved are in the order of a 2 , aa s and a 2 , respectively. 
Consequently the corresponding events must have quite different weights. In the 
total cross section, however, the order of a is common for these three channels. In 
other words the probabilities that a single real photon resolves or not is 0(a) : 
0(1). Hence we can rearrange a by removing it from the hard cross section to the 
probability of the photon branching. Thus we have the ratio 0(a) : 0(a) to branch 
or not for a photon and 0(1) : 0(a s ) : 0(a 2 ) for the hard cross sections, which are 
comparable in magnitude. This rearrangement then allows us to generate events 
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with good efficiency. We use da jdp\ = n/p^ as the reference cross section. If any 
hard cross section is equal to this, the event is generated with weight = 1. 

We like to point out also that the argument of a s in the hard cross section is p^, 
not s in our model(in contrast to the case of QCD branching where the virtuality K 2 
enters into the argument). This is because there appear some t-channel diagrams in 
the hard scattering, and when p\ is smaller than s, the cross section is dominated 
by the t-channel contribution. The running effect comes from the virtual corrections 
in the t-channel propagator and its scale is developed to p\. 

When the final states of the hard scattering is fixed, we make timelike evolutions 
from Q 2 for timelike partons. Then we fix the energy of the initial photons in 
the electron and the positron, and thus we can determine the four-momenta of all 
partons. 

Finally a comment is given about the higher order contributions. In the literature 
it has been emphasized that they are important for quantitative analysis |]TT], |12 



It has been known, however, from the study of e + e~ annihilation, that the most 
important contributions in the next-to-leading logarithmic terms can be included in 
the LL shower by using the transverse momentum squared p\ as the argument of 
the QCD coupling a s and imposing the angular ordering]?!. This is also valid for 
the photon evolution we are considering. These effects are easily implemented into 
our model. 



Section 4. Results 

We present results on parton distributions. These distributions are usually ob- 
tained by integrating the AP equation numerically. The parton shower provides us 
an another method to get them. Though this equivalence has been well established, 
it will be instructive to understand the algorithm, by comparing the results of the 
parton shower with those obtained by the conventional method. Many kinds of 
distributions JT^] have been proposed in the literature. First we use those given in 
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ref. |nj for the LL approximation, because their distributions are parametrized in 
simple functions. In order to make precise comparison, we calculate parton distribu- 
tions by the shower model with a fixed cutoff e to regularize the infrared singularity. 
Results are shown in Fig.2, where one can see good agreement which shows that 
the shower algorithm yields the same distribution as obtained by usual methods to 
solve the AP equation directly. 

Next we discuss the higher order effects. As mentioned in section 3, the most 
important effect beyond the leading order can be easily included into the LL model 
by using the transverse momentum squared p\ as the argument of the running 
coupling. To show it, we calculate the structure function, F 2 (x, Q 2 ) by the parton 
shower with the transverse momentum p\ and the virtual mass K 2 as the arguments. 
As the full next-to-leading logarithm(NLL) solution we take the distributions given 
in [O. Here one should note that this comparison is not trivial. One reason is that 



the parton distributions in the NLL order are scheme-dependent and the definition 
of the initial distributions is not unique. Another reason is that in ref. [12] the 
evolution started at too small Ql so that contributions beyond the NLL order cannot 
be ignored. Therefore for the comparison we start the shower evolution at Q 2 = 
4GeV 2 using F 2 (x, 4GeV 2 ) as the initial parton distribution. From Fig.3 we find 
that our shower model with the running coupling of p\ can reproduce the NLL 
results qualitatively. This agreement is more remarkable for higher Q 2 . However 
considering the fact that we are able to get a qualitative agreement between the 
NLL calculation and the LL shower together with the fact that there still remains 
non negligible experimenatal errors in the present data, there seems no need for the 
NLL shower model for the moment, though it is in principle possible 0. 

Next we will discuss jet production in two quasi-real photon processes in e + e~ 
scattering. In order to check the program, we compare the results of our model with 
the inclusive QCD calculations obtained by integrating the parton distributions 
fi(x,PT) with the hard cross sections , where i,j denote photon, quarks and 
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gluon. 

da 



dp\ 



( dziF l/e {z u E) ( dz 2 F l/e (z 2 ,E) 

JO « 



/•i /•! da-(s) 

dxxfifa, Pt) / dx 2 fj(x 2 , Pt)S(zi z 2 x x x 2 s - s) ? J 2 (14) 

■ ■ JQ J Ci'Prp 



I. J 



It should be noted that it is not trivial for our model to get the cross section, 
because the parton distributions are obtained as a result of the generated events |L5 



not being assumed at the beginning of the calculation. Figure 4, where the inte- 
grated cross sections, J 2 dp\da jdp\ are given, shows consistency of our model. 
Here to make a detailed comparison, we present the results for hard scattering 
of quark-photon and gluon-photon. The contributions from other processes are 
small (parton-parton scattering) in TRISTAN energy range, or are trivial (photon- 
photon scattering). 

The parton shower model has an advantage, when one wants to make topological 
analyses of the final particles. As a typical example of these kinds of analyses, we 
take the thrust distribution. Since the system of final partons is boosted along z- 
direction, we analyze events using transverse momentum. Similarly to the case in 
e + e~ annihilation, the thrust distribution changes as the energy increases(see Fig. 5). 
Note again that in the present analysis we make no hadronization and the results are 
not fully meaningful as a physical prediction. The hadronization should be included 
for quantitative analysis of actual data. 



Section 5. Summary and discussions 

In this paper we proposed an algorithm of parton shower model for two-photon 
processes in e + e~ scattering and constructed a Monte Carlo model according to it. 
We formulated the algorithm using AP equations and the kernel function method || 
for the QCD evolution. In the construction of the model, we adopted the forward 
spacelike evolution for the energy distributions of partons xq(x, Q 2 ) by taking ac- 
count of the breaking of energy conservation due to the inhomogeneous term in the 
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AP equations for the photon. 

Using this Monte Carlo model, we have shown the parton distributions of the 
photon and some results on jet production in the two-photon process. These results 
clearly demonstrated that the model is consistent with perturbative QCD. 

Our model will be completed as an event generator if combined with some 
hadronization models. Then it can be used for reliable estimation of the detec- 
tor efficiency in experimental analysis. On the other hand the model will be useful 
for theoretical investigations, as the model is rigorously based on perturbative QCD. 
For example, it is easy to change the input parton distributions at Ql in our model. 
Thus we can study in detail how the change of the input distributions affects the 
final results. Also interesting is the study on charm distributions. However, the 
value of Qq c , where the photon starts to resolve into charm quark pair, is not clear. 
Naively one expects 4m^, but kinematically it should be Am 2 c / xj,, where Ob}) IS db 
Bjorkcn variable. The parton shower model presented in this work would give us 
a quantitative answer to the question whether this kinematical change modifies the 
parton distributions or not. Extensive studies on these points will be presented in 
a forthcoming paper. 
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Figure Captions 

Fig.l The schematics of deep inelastic scattering of photon(a) and two quasi- real 
photons process(b) in e + e~ scattering. Here 7 denotes a quasi-real photon, while 7* 
denotes the virtual one. 

Fig.2 Parton distributions xf(x, Q 2 )/a. The solid and dashed curves are calcu- 
lated by the parametrized form given in ref. [14] for LL approximation. The circles 
and diamonds show results of our shower model with a constant cutoff, which should 
be very small compared to the observed (1 — x) value. Here it is 10~ 4 . The evolu- 
tion starts from Q 2 = 0.25GeV 2 by taking common initial distributions, (a) w-quark 
distributions, (b) -u-quark and gluon distributions. 

Fig.3 The structure function F 2 (x,Q 2 )/a. The curves are calculated by the 
authors of ref. [12] The diamonds are results by the shower model with p 2 ^, while 
the plus symbols represent that with the virtual mass squared. These distributions 
are normalized at Q 2 = 4GeV 2 . The full circles are the experimental data[l]. Q 2 = 
73GeV 2 (a), 1000GeV 2 (b). 

Fig.4 Cross sections of jet production integrated over the transverse momentum 
in the two quasi-real photon process in e + e~ scattering. The solid and dashed curves 
are calculated with the parton distributions of ref. [14]. Circles and diamonds are 
results of our parton generator. In the solid curve and circles the hard scattering 
contains only the photon-quark scattering, while in the dashed curve and diamonds 
it contains only photon-gluon scattering. The c.m.s energy is 58GeV. 

Fig. 5 Thrust distributions of the final partons in the two quasi- real photon 
process in e + e~ scattering. This is calculated only by transverse momenta of 
the final partons. The histogram, circles and diamonds are results with Pr min = 
20GeV 2 , 50GeV 2 and lOOGeV 2 , respectively. The c.m.s energy is 58GeV. 
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